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T.  INTRODUCTION 

When  a  chemical  reaction  takes  place  old  chemical  bonds  are  broken  and 
new  ones  formed.  The  energy  associated  with  these  bonds  either  appears  in 
the  form  of  heat  or  results  in  the  absorption  of  heat  from  the  surroundings. 

A  rapid  exothermic  reaction  in  which  some  of  the  reactants  are  oxidized  is 
*  'town  as  a  combustion  reaction.  The  zone  in  which  combustion  reactions  take 
place  is  known  as  the  reaction  zone  and  as  it  moves  relative  to  an  unburned 
Fixture  it  is  known  as  a  combustion  wave.  . 

“  *  r  1  ‘  * ~r  ~  it 
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A  detonation  wave  is  a  special  case  of  the  combustion  wave  in  that  it 
r’oves  with  a  supersonic  velocity  relative  to  the  unburned  mixture  and  it  is 
sustained  by  the  energy  liberated  in  the  chemical  reaction.  It  is  composed 
uf  a  shock  wave,  or  discontinuity,  which ’propagates  into  an  explosive  material, 
expressing  and  heating  it,  thus  triggering  chemical  reaction.  A  balance  is 
achieved  such  that  the  chemcial'  reaction  then  supports  the  shock. 

Detonation  waves  have  been  studied  extensively  from  both  experimental 
;.n !  theoretical  points  of  view  for  several  years.  Considerable  work  tnnr 
been  done  on  the  detonation  of  explosives  that  exist  in  the  gaseous,  liouid, 
conJensed  solid  state.  On  the  other  hand,,  a  ‘  1  imi ted  amount  of  attention 
-a;  been  given  to  the  study  of  a  cloud  of  powdered  explosive  dispersed  in 

i  *  r  m 
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II.  THEORETICAL  DESCRIPTION 

After  initiation,  the  detonation  will  reach  a  steady  state  in  which  it 
propagates  into  an  unreacted  mixture  at  a  constant  supersonic  velocity.  In 
order  to  predict  the  pressures  generated  by  such  a  wave,  among  other  things, 
it  is  necessary  to  formulate  the  phenomenon  analytically  as  was  first  done 
!.y  D.  Chapman  and  E.  Jouget  around  1900.  Results  of  their  work,  known  as 
C.1  theory,  are  still  utilized  extensively  today. 

•  ,v  •  *>■' .  *  ■  . 

•  A.  Chapman  -  Jouget  Theory 

i 

T!.e  assumptions  made  in  order  to  formulate  the  CJ  model  of  a  detonation 
•»-we  jre  as  follows:  '  .  . 

1.  The  flow  is  one-dimensional  and  laminar. 

*  •  i 

2.  The  detonation  wave  is  pictured  os  a  jump  discontinuity,  a 

» 

'hock  across  which  the  basic  conservation  laws  apply.  The  chemical  reaction 
't'larded  as  instantaneous  so  that  there  is  essentially  no  reaction  zone, 
only  a  region  in  whic.i  the  unburned  mixture  exists  in  thermodynamic  equilibrium 
and  a  . region  in  which  the  gaseous  products  of  reaction  exist  (also  in  therrro- 
o .n.unic  equilibrium),  the  regions  separated  by  the  shock. 

3.  The  motion  of  the  detonation  as— i t  moves  into  the  explosive  is 
vrependent  of  time.  (Steady) 


0 


B.  Basic  Conservation  Laws 


If  we  consider  a  detonation  wave  moving  with  a  velocity  D  into  an  explosive 
initially  at  rest  then  the  basic  conservation  laws,  when  applied  to  such  a 
system,  result  in  the  following  equations. 


Conservation  of  Mass:  pQ0  =  p(D-U)  (1) 

where  =  density  of  the  unreacted  explosive  (initial  state) 

p  3  density  of  the  gaseous  product  mixture  (final  state) 

U  =  velocity  of  the  gaseous  products 
*  •  ,  •  •  ■ 

2  ? 

Conservation  of  Momentum:  pQ  +  pQ0  =  p  +  p(D-U)  (?) 

where  pQ  =  pressure  of  initial  state 
p  3  pressure  of  final  state 

'  "nservation  of  Energy:  E  +  i  (D-U)2  =  EQ  +  +  1  qc>  (3) 

where  E  =  specific  internal  energy 


:im  t  ion  of  Sta  te  :  E  =  E  ( p.p  ) 

Fcr  an  ideal  gas 

E  -  E  3  C  T  -  C  T-  -  q  =  -  (  ^.)  -  1°  ( 

o  v  vo  o  ^  v  y-1 '  ~  v 

o 

where  Cy  =  constant  volume  heat  capacity 
q  3  heat  of  reaction 
r  -  ratio  of  specific  heats 


& 


-  o 


(•1) 


de  have  therefore  a  system  of  four  equations  in  the  five  unknowns 
.  ,  o,  end  q.  Further  assumptions  must  therefore  bo  made  in  order  to  a 
state  of  the  reaction  products  following  the  dc ton j* :on . 


»  > 


C.  Determination  of  the  CJ  Detonation  Velocity 

It*  U  is  eliminated  from  Eqs.  (1)  and  (2),  one  obtains  an  equation  that 
r.«>  be  plotted  as  a  line  in  the  p-v  (v-l/p)  plane  and  is  known  as  a  Rayleigh 
Inc*.  This  lino  shows  the  variation  of  product  pressure  and  specific  volume 
tor  particular  detonation  velocities.* '  (Fig.  la) 

If  U  and  D  are  eliminated  from  Eqn.  (3)  via  Eqns  (1)  and  (2),  one  obtains 
.in  equation  that  may  be  plotted  in  the  p-v  plane  and  is  known  as  a  Hugoniot 
I'.n-vi*.  This  gives  the  volume  as  a  function  of  internal  energy,  E  (p,v),  which 
is  given  as  an  equation  of  state'.  (Fig.  lb) 

I 

The  state  (p,v)  of  the  reacted  products  may  than  be  determined,  for  a 
•  i ' r‘ n  detonation  velocity ,  as  the  intersection  of  the  Rayleigh  line  an.) 

H.rj'liot  curve.  (Fig.  1c) 

t  • 

The  expression  for  the  Rayleigh  line  contains  the  detonation  velocity  as 

a  ’’arametfir  so  the  entire  family  of  Rayleigh  lines  should  be  considered,  three 

of  which  are  shown  in  Fig.  lc. 

» 

For  a  sufficiently  small  detonation  velocity,  say  D.,,  the  two  curves  do 
in.;,  intersect  so  there  is  no  solution.  For  ajarge  value,  ,  there  are  two 
solutions,  denoted  S  (strong)  3nd  VI  (weak).  And  for  a  particular  value  of  0, 

D  ,  cnere  is  one  unique  solution.  The  Chapman  -  Jouguet  hypothesis  is  that  the 
:*a!y  detonation  velocity  is  indeed  this  value,  i.e.,  the  minimum  velocity 
'  -.n*,  is  tent  with  the  conservation  equations.  Other  solutions,  for  D>Dj  are 


'Vi 


rejected  because  of  physical  reasonings. 

0.  Solution  ot  Governing  Equations  Under  the  CJ  Hypothesis 

The  CJ  hypothesis  may  he  used  together  with  one  further  assumption  in 
order  to  obtain  a  very  simple  system  of  governing  equations.  This  assumption 
is  essentially  that  the  detonation  wave  is  a  “strong"  reaction  or  disturbance. 
In  other  words,  one  might  also  say  that  the  Mach  number  of  the  detonation  wave 
is  large  or  that  the  pressure  change  resulting  from  the  passage  of  the 
detonation  wave  Is  large.  If  such  is  the  case  then  the  system  of  governing 
equations  may  be  reduced  to 

1 

.  .  V  7  2  <5> 

o  •  e0  (  *  p  “  ~t*T  0  ' 

If,  therefore,  the  density  of  the  Initial  explosive,  the  heat  of  reaction,  and 
tie  product  ratio  of  specific  heats  are  known,  the  CJ  pressure,  density,  and 
detonation  velocity  may  be  determined.  ;  It  Is,  however,  a  rather  formidable 
task  to  compute  the  equilibrium  composition  of  the  combustion  products  (which 
r;  necessary  to  obtain  y  for  example),  and  this  r.  almost  always  handled  by 
means  of  established  computer  codes. 


■T- 


III.  PURPOSE  OF  RESEARCH 

Two  groups  have  been  involved  in  theoretical /experimental  work  on 
the  problem  of  the  detonation  of  dispersed  high  explosives:  The  Gas 
Dynamics  Laboratory  at  the  University  of  Michigan  and  the  Illinois 
Institute  of  Technology  Research  Institute  (IITRI)  Both  groups 

utilize  the  CJ  detonation  model  to  theoretically  predict  the  pressures 
and  other  flow  variables.  The  governing  equations,  Eqs.  (5),  indicate 
that  it  is  the  composition  of  the  products  of  combustion  that  influence 
the  CJ  pressure,  detonation  velocity,  etc.  It  is  therefore  necessary  to 
determine  the  equilibrium  state  of  the  combustion  products.  The 
Michigan  group  utilizes  the  Gordon -McBride  computer  program  to  do 

so,  whereas  IITRI  uses  their  own  EQUIL  code 

It  is  the  purpose  of  this  report  to  discuss  three  established 
computer  codes  that  are  used  in  the  calculation  of  equilibrium 
composition  of  combustion  products,  namely: 

A.  Gordon  and  McBride's  program  as  described  in  NASA  SP-273, 

B.  EQUIL  program,  written  by  IITRI,  and 

C.  TIGER,  written  by  Ballistic  Research  Laboratories 


IV.  DESCRIPTION  OF  COMPUTER  CODES 

A.  Gordon  and  Mc8ride 

The  Gordon -McBride  code  is  a  well -accepted  and  readily  available 
code  that  can  be  utilized  to  calculate  equilibrium  compositions  for  CJ 
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detonations.  The  program  combines  the  thermodynamic  equilibrium 
calculations  and  the  CJ  Equations  (Eqs.  (5))  to  completely  define  the 
detonation  wave.  The  condition  for  equilibrium  may  be  stated  in  terms 
of  any  of  several  thermodynamic  functions  such  as  the  free  energy  or 
entropy.  This  treatment  iteratively  minimizes  the  free  energy  to 
calculate  the  equilibrium  composition  and  properties  of  perfect  gas 
mi xtures. 

There  are  three  steps  in  the  calculation  of  the  CJ  parameters. 

First  an  initial  estimate  is  made  of  the  detonation  pressure  and 
temperature.  Then  an  improved  estimate  is  obtained  via  a  recursion 
formula  and  finally  the  correct  values  evaluated  by  a  Newton-Raphson 
iteration  procedure. 

The  input  specifies  the  fuel  and  oxidizer  composition,  enthalpy, 
density,  and  the  equivalence  ratio.  The  output  includes  properties  of 
the  burned  and  unburned  gases,  detonation  velocity,  detaonation  ratios 
P/Po,  T/To,  fVp0»  and  the  final  equilibrium  composition. 

B.  EQUIL  (8). 

The  computer  program  EQUIL  was  written  by  Dr.  R.  Snow  for  IITRI  and 

is  "proprietary  material"  and  not  obtainable,  hence  thepresent 

discussion  is  based  solely  on  information  supplied  by  that  contractor. 

* 

It  can  be  shown  that  if  Pv  =  pressure  reached  in  a  constant  volume 
explosion  then  PCJ.  =  2Py.  The  IITRI  EQUIL  computer  program  is  based  on 
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the  minimization  of  Gibb's  free  energy,  which  should  be  done  under  the 
constraint  of  constant  temperature.  By  using  various  heat  balances, 
however,  the  program  computes  the  pressure,  not  at  constant  temperature 
but  at  constant  volume,  and  then  relates  this  value  to  the  CJ  pressure 
through  the  above  relationship. 

C.  TIGER 

TIGER  is  a  highly  complex  computer  code  designed  to  calculate 
detonation  properties  under  a  wide  variety  of  thermodynamic  conditions, 
for  example,  utilizing  different  equations  of  state.  The  four  volume 
documentation  listed  in  the  bibliography  includes  theoretical 
discussions  and  will  be  referenced  here  rather  than  discussed  in  detail. 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 

All  three  computer  programs  are  dependent  upon  thermodynamic  data 
as  input  for  their  calculations.  The  acquisition  of  acceptable  data  for 
explosive  compounds  is  a  major  task  in  itself  and  numerous  sources  are 
presently  available  (for  example  see  Refs.  (1),  (3),  (7)).  Conflicting 
values  are  present  in  almost  all  sets,  and  there  are  numerous  arguments 
in  the  current  literature  as  to  why  one  version  is  superior  to 
another.  The  JANAF  tables  ^  have  long  served  as  a  standard  and, 
outside  of  that  source,  it  is  difficult  to  obtain  values  for  S°  for 
example. 

Using  the  JANAF  tables,  one  can  run  the  Gordon-McRride  program  for 
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various  fuel/oxidants  and  initial  conditions.  A  typical  output  is 
attached.  This  program  is  well -documented,  easily  run,  and 
although  limited  to  perfect  gases,  gives  feasible  results  that  could 
provide  a  starting  point,  should  further  refinements  become  necessary. 

The  EQUIL  program,  since  it  is  proprietary  in  nature,  cannot  be 
obtained  and  one  must  rely  on  IITRI  and  the  corresponding  interpretation 
of  results.  This  has  been  done  with  confidence  in  the  past  and  it 
is  felt  that  the  two  programs  will  give  results  that,  although  not  in 
exact  agreement  by  any  means,  differ  by  an  acceptable  amount,  if 
corresponding  thermodynamic  data  are  used.  It,  also,  is  limited  to  the 
case  of  perfect  gases. 

TIGER,  written  in  FORTRAN  IV,  is  an  immense  program  that  is  neither 
user-friendly  nor  system-f riendly .  Large  amounts  of  man-hours  and 
computer  funds  were  expended  and  no  acceptable  results  obtained.  It 
supposedly  has  the  capability  of  providing  various  equations  of  state 
and  many  other  refinements,  yet  is  a  code  that  the  somewhat -casual  user 
should  approach  with  extreme  caution. 

In  conclusion,  as  one  seeks  to  numerically  determine  the 
equilibrium  composition  of  combustion  products,  one  must  first  have 
acceptable  thermodynamic  data,  regardless  of  which  code  is  used,  since 
each  and  every  code  is  highly  dependent  upon  such  data.  Gordon  and 
McBride  offer  a  readily  available,  user-friendly  code  that  produces 
results  which,  although  limited  in  many  ways,  provide  an  excellent  first 
plateau.  For  many  applications  no  further  refinements  would  be 
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necessary.  The  EQUIL  program  is  greatly  limited  because  of  its 
availability  although  results  obtained  from  that  program  have  proved 
satisfactory  when  compared  with  experiment  .  The  TIGER  program 
offers  the  highest  degree  of  flexibility  and  power,  but  one  must  be 
intimately  associated  with  the  problem  under  consideration  and  the  code 
itself  in  order  to  avail  oneself  of  such  power. 
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